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ABSTRACT 

We spectroscopically survey the galaxy cluster XMM-LSS J02182-05102 (hereafter IRC 0218) using 
LRIS (optical) and MOSFIRE (near-infrared) on Keck I as part of the ZFIRE survey. IRC 0218 has a 
narrow redshift range of 1.612 <Zspec< 1.635 defined by 33 members of which 20 are at Rproj< 1 Mpc. 

The cluster redshift and velocity dispersion are Zci = 1.6233 ± 0.0003 and tJci = 254 ± 50 km s“*. 

We reach NIR line sensitivities of ^ 0.3 xl0“*^ erg s“*^ cm“^ that, combined with multi-wavelength 
photometry, provide extinction-corrected Ha star formation rates (SFR), gas phase metallicities from 
[Nil]/Ha, and stellar masses. We measure an integrated Ha SFR of ~ 325 Mq yr“* (26 members; 

Rproj< 2 Mpc) and show that the elevated star formation in the cluster core (Rproj< 0.25 Mpc) is 
driven by the concentration of star-forming members, but the average SFR per Ha-detected galaxy 
is half that of members at Rproj"^ 1 Mpc. However, we do not detect any environmental imprint 
when comparing attenuation and gas phase metallicities: the cluster galaxies show similar trends with 
M* as to the field, e.g. more massive galaxies have larger stellar attenuation. IRC 0218’s gas phase 
metallicity-M* relation (MZR) is offset to lower metallicities relative to z ^ 0 and has a slope of 
0.13 ± 0.10. Comparing the MZR in IRC 0218 to the COSMOS cluster at z = 2.1 shows no evolution 
(At ^ 1 Gyr): the MZR for both galaxy clusters are remarkably consistent with each other and 
virtually identical to several field surveys at z 2. 

Subject headings: galaxies: clusters: individual (XMM-LSS J02182-05102) - galaxies: evolution - 
galaxies: abundances - galaxies: fundamental parameters 


1. INTRODUCTION 

The advent of sensitive infrared imaging surveys has 
provided a powe rful method for identifying galaxy clus¬ 
ters a t z > 1.5 (lEisenhardt et al.l 120081: l^povich et al.l 
120101 : iGonzalez et al.l l2012t iSoitler et al.l I2012D . We 
are now able to track how clusters build up their 
galaxy populations over ~ 80% of cosmic time 
(O < z < 2.1). Of particular interest to cur¬ 
rent galaxy formation models are how scaling relatio ns 
evolve and depend on environment (iDave et al.ll2011bllal : 

* http://zfire.swinburne.edu.au 

^ George P. and Cynthia W. Mitchell Institute for Funda¬ 
mental Physics and Astronomy, Department of Physics & As¬ 
tronomy, Texas A&M University, College Station, TX 77843; 
vy@physics.tamu.edu 

® Swinburne University of Technology, Hawthorn, VIC 3122, 
Australia 

Research School of Astronomy and Astrophysics, The Aus¬ 
tralian National University, Cotter Road, Weston Creek, ACT 
2611, Australia 

® Astronomy Department, Yale University, New Haven, CT 
06511, USA 

® Mitchell Postdoctoral Fellow 

^ Department of Physics and Astronomy, The University of 
Kansas, Malott Room 1082, 1251 Wescoe Hall Drive, Lawrence, 
KS 66045, USA 

® Department of Physics and Astronomy, University College 
London, Gower Street, London WCIE 6BT, UK 

® Department of Physics and Astronomy, Faculty of Science 
and Engineering, Macquarie University, Sydney, NSW 2109, 
Australia 

Australian Astronomical Observatory, PO Box 915, North 
Ryde, NSW 1670, Australia 

Leiden Observatory, Leiden University, P.O. Box 9513, NL- 
2300 RA Leiden, The Netherlands 


Genel et all 12014 iTonnesen fc Ceiil 12014 iSchave et al.l 
20151 iTavlor &: Kobavashill2015n . For example, the ac¬ 
tivity observed in clusters at z > 1 combined with the 
old stellar ages of cluster galaxies at z 0 indicate that 
1 < z < 2 is an important epoc h for understanding 
how star formation i s quenched fe.g. lRudnick et ani2012l : 
iBrodwin et al.l[20T^ . 

Essential to measuring scaling relations at z ^ 2 is 
separating first cluster galaxies from the field. IR imag¬ 
ing is effective at identifying galaxy overdensities at 
z > 1.5, but broadband photometric redshifts (zphot) are 
coarse and field contamination is a concern, especially 
beyond the cluster cor43. This is also true for stud¬ 
ies using narrow-band imaging to identify Ha-em itting 
members (|Finn et al.lI2005t iKovama et aLll2013aF) . Nei¬ 
ther method can measure cluster kinematics, and both 
are susceptible to contamination by broad-line Active 
Galactic Nuclei (AGN). While Hubble Space Telescope 
WFC3 grism observations can measure continuum fea¬ 
tures, e.g. the Balmer break, a s well as strong emis¬ 
sion lines (IBrammer eF "all 120121 IZeimann et al.l 120131 
INewman et al.ll20l4IGobat et al.ll2013ll . the spectral res¬ 
olution is too low for measuring cluster kinematics and 
spectral lines such as Ha and [Nil] are blended. The 
WFC3 footprint also is limited to the cluster core. 

Due to these limitations, currently there are only four 
(proto) galaxy clusters at z > 1.5 that have enough 
members (> 15) for a kinematic analysis (|Bavliss et al.l 

Although see results using medium-band NIR filters that im- 
prove 2phot by a factor of 5 compared to broadband measure¬ 
ments JWhitaker et al.]|2011l : [5pitler et al.ll201^ 
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12014 lYuan et all 12014 iShimakawa et alJ 12014 . The 
large range in measured cluster velocity dispersions (550- 
1500 km s“^) and substantial spatial substructure con¬ 
firm that these are dynamically young systems. Because 
the z ~ 1.5 — 2 clusters are still assembling, we can 
better disentangle evolution driven by environment ver- 


sus galaxy mass ( 

Pens- et alJ I201GI: IWetzel et al.l 120121: 

Ouadri et al.l 12012 

: Muzzin et al. 120121: IPanovich et al.l 

20121: iBassett et al.ll2013ll. For example, the discovery of 


intense star formation in the cluster cores (iHilton et al.l 
l2010t iTran et ^120101 : IZeimann et al.ll20Tll suggests an 
active interplay between enrichment of gas in the Inter- 
Stellar Medium (ISM) and the Intra-Cluster Medium 
(ICM). 

The primary challenge to studying cluster galaxy 
populations at z > 1.5 is the need for near-lR 
spectroscopy. With the installation of efficient multi- 
obiect near-lR spect rographs such as Subaru/FMOS 
(IKiinura et al.l 1201011 . VLT /K MOS (ISharples et al.l 
12014 . and Keck/MOSFIRE (IMcLea.n et al.l 1201 2n . 
we can now obtain key rest-frame features from 
[Oii]A 3727,3729A to [Sii]A6717,673lA for galaxy 
clusters at z > 1.5. These spectral features have been 
used to establish empirical scaling relations at z ^ 0 
that we can now test at z ^ 2, i.e. can we detect 
an environmental dependence in addition to redshift 
evolution? 

The increasing activity in galaxy clusters with red- 
shift raises the question of how star formation rates 
(SFR), stellar masses (M*), attenuation, and gas phase 
metallicities are related. Star-forming cluster galaxies 
at z > 1 define a SFR-M,^ relation similar to the field 
(“main sequence”; iDaddi et al.ll2007llNoeske et al.ll2007l : 

iKovama et al.ll2013b[ l. but studies suggest there are envi¬ 

ronmental differences, e.g. suppressed st ar formation in 
low-m ass members and less attenuation (jZeimann et al.l 
1201311 . If cluster galaxies follow a different relation be¬ 
tween specific SFR and stellar mass (SSFR-M*), this 
could mean that the ioniz ation conditions are also dif¬ 
ferent (jKewlev et al.ll201^ . 

Recent studies of gas phase metallicity versus M* 
(MZR) find that galaxie s at z ~ 2 are unifornily offset 
from the local relation (Erb et all l2006at iSteidel et al.l 
12014 iSanders et al.l[2M5 1. However, while the MZR in 
the local universe shows an environmental dependence, 

(iCooper et al.l[2008l : lEllison et al.l[2009l : lPeng fc Maiolinol 

12014 . it is unclear if this is also the case at z ^ 
2: several studies using stacked spectra find the clus¬ 
ter galaxies to be more meta l -rich than their field 
counterparts llKulas et al.l 12013: IValentino et al.l 12014 
IShimakawa et al.l l2015h . but Kacprzak et al.l (l2015fl 
study ~ 50 cluster galaxies at z = 2.1 and find no differ¬ 
ence from the field. 

To determine if empirical relations vary with environ¬ 
ment at z ^ 2, the ZFIRE survey combi nes NIR spec- 
troscopy obtained with Keck/MOSFIRE (iMcLean et al.l 
1201211 with deep multi-wavelength imaging to study clus¬ 
ter galaxies at z ^ 2. We compare to held measure- 
ments obtained by c omplementary su rveys such as KBS S 
(|Steidel et al.ll20l4 and MOSDEF (|Kriek et al.ll2015ll . 
Initial ZFIRE results on the COSMOS cluster at z = 2.1 
(|Spitler et al.l 1201211 include spectroscopically conhrm- 
ing > 50 members and measuring the cluster’s veloc¬ 


ity dispersion (|Yuan et al.l 12014 , mea suring the clus¬ 
ter’s gas phase metallicity-M* relation (|Kacprzak et al.l 
1201511 . and investigating the ionization properties of the 
Inter-Stellar Medium (ISM) of individual cluster galaxies 
(|Kewlev et al.l[2M^ . 

As part of our ZEIRE survey, we obtain optical and 
NIR spectroscopy of XMM-LSS J02182-05102 (here¬ 
after IRC 02181^ at Zni = 1. 6233. Originally re¬ 
ported by iPapovich et al.l (I2010I1 . IRC 0218 was one 
of the hrst galaxy clusters to show an increasing frac- 
tion of star form ation with increasing local density 
(|Tran et al.ll2010ll . Estimates of IRC 0218’s total star 
formation rate as determined by IR observations are > 
1000 Mq yr“^ (|Popesso et al.l[20Tl : ISantos et al.ll20l4 . 
Because IRC 0218 is still in the proce ss of building its 
galaxy population (|Rudnick et al.l[20T^ , the cluster is an 
ideal target for tracking how star formation, gas phase 
metallicities, and dust vary with stellar mass as a func¬ 
tion of environment at z > 1.5. 

In our analysis, we center IRC 0218 on the Bright¬ 
est Cluster Galaxy (BCG) at (a, S)j 2 ooo = (02 : 18 : 
21.5, —5 : 10 : 19.9). The BCG is also currently the most 
distant galaxy-galaxy lens with a total mass withi n the 
Einstein radius of 1.8 x 10^^ Mq (iWong et al.ll20l4 . We 
use a Chabrier Initial Mass Function and AB magnitudes 
throughout our analysis. We assume = 0.7, Ha=0.3, 
and Hq = 70 km s“^ Mpc“^. At z = 1.62, the angular 
scale is 1" = 8.47 kpc. 

2. OBSERVATIONS AND DATA 
2.1. Target Catalog 

Spectroscopic targe ts were selected from the 
iWilliams et al.l (1200911 catalog of the Ultra-Deep 
Survey (UDS) taken as part of the UKIRT Infrared 
Deep Sky Surv ey (UKIDSS), a pu blic near-infrared 
imaging survey (|Lawrence et al.l [200711 . The AT-selected 
catalogs reach 5cr-limiting magnitudes in a 1.75" diame¬ 
ter of Hab < 27.7, i?AB < 27.1, zab < 26.8, zab < 25.5, 
Ja b < 23.9, and A I ar < 23.6 mag. We refer the re ader 
to IWilliams et al.l (1200911 and iQuadri et al.l (1201211 for 
more details on the photometric and Zphot catalogs used 
in this work. 

2.2. Optical Spectroscopy: Keck/LRIS 

Usin g! the Low Re solution Imaging Spectrometer 
fLRIS: ^ke et al.lll995L 5.5' x 8' field of view; 0.135" 
per pixel) on Keck, we carried out a spectroscopic sur¬ 
vey centered on the cluster on 19 & 20 October 2012 
(NASA/Keck Program ID 48/2012B). The primary LRIS 
targets were selected to be brighter than zab= 21 
mag and the secondary targets (mask fillers) between 
21 <zab< 24 mag. No morphological selection was ap¬ 
plied. Of the primary targets, higher priority was also 
given to objects i dentified to be ca ndidate star-forming 
cluster galaxies (ITran et al.l 1201011 . candidate Lyman- 
Break Galaxies at Zphot> 1.35, a nd [Oii]-emitters id en- 
tified from narrow-band imaging (|Tadaki et akllMT^ . 

We use the 600/4000 grism for the blue side of LRIS 
(0.38 < A < 0.58 ^m) and the 600/10000 grating for the 
red side (0.70 < A < 1.0 ^m); with 1" slit widths, the 

This galaxy clust er is also refe r red to as CLG0218.3-0510 by 
I'llan et al.l II2010I1 and ISantos et al.l II2014I 1. 
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corresponding resolution is 4.0A and 4.7A respectively. 
Observing conditions were excellent with median seeing 
of about 0.6". A subset of high priority targets, e.g. the 
brightest cluster galaxies, were targeted in three of the 
four observed masks. For masks 1, 2, and 3 which in¬ 
cluded common high priority targets, we obtained 9 x 20 
minute exposures. For mask 4 which included a more 
general redshift selection, we obtained 5 x 20 minute ex¬ 
posures. The four LRIS masks have a total of 136 targets 
including repeats. For all the exposures, we dithered by 
±1" along the slit to minimize sky residuals. Standard 
calibrations including flat-fields and arc lamps (ArCd- 
HgNeZn) for each mask also were taken. 

To reduce the spectra, we follow pRan et al.l (1200711 and 
use IRAF F^ routines with cust om software provided by 
D. Kelsoii ilKelsonl [19981120(3^ . To summarize, we first 
convert the multi-extension files into single images for the 
blue and red sides; we reduce the blue and red sides of 
LRIS separately. Cosmic rays are particularly problem¬ 
atic on the red side and we remove as many as possible 
using crutil in IRAF; the blue side is minimally affected 
by cosmic rays. We flat-field the science frames, deter¬ 
mine the vacuum wavelength solution with the arclamp, 
remove the sky lines, and rectify the 2D spectra. We me¬ 
dian the rectified science frames for each mask to obtain 
a single combined image. 

The 2D spectra first were visually inspected to identify 
all potential sources including emission-line only objects. 
The ID spectra were extracted by summing the five rows 
centered on the source (0.135" per pixel) of each 2D spec¬ 
trum; the wavelength coverage for the extracted spectra 
is 3800 —5800A (blue side) and 7000 —lOOOOA (red side). 
For galaxies at z ~ 1.6, this corresponds to rest-frame 
UV where the continuum is detected for only a hand¬ 
ful of systems. To mea sure redshifts, we use xcsao in 
IRAF (jKurtz et al.lll99^ with various templates for star¬ 
forming, quiescent, post-starburst, Lyman-break, and 
Lya galaxies; we refer the reader to iTran et al.l (|2005bf l 
for details. 

2.3. Near-Infrared Speetroscopy: Keck/MOSFIRE 

At Zci = 1.6233, most of the rest-frame optical spec¬ 
tral features such as H/3 and Ha have shifted to the 
near-lR. To measure these features and provide contin¬ 
uous spectral cov erage of the cluster galaxies, we used 
Keck/MOSFIRE ([McLean et al.llMll on 23 & 24 De¬ 
cember 2013 and on 10-13 February 2014. Primary tar¬ 
gets were galaxies spectroscopically confirmed with the 
LRIS observations to be at z > 1.5 (see Fig. [1]). We also 
included objects from the 3D-HST survey with grism red- 
shifts z > 1.5 (Momcheva et ah, in prep). Secondary tar¬ 
gets were selected to have Zphot> 1-4 using the Williams 
catalog, and the lowest priority fillers had Zphot< 1-4 and 
fAB< 22 mag. 

We obtained YJH spectroscopy corresponding to wave¬ 
length ranges of 0.97 — 1.12/rm, 1.15 — 1.35/rm, and 
1.46 — 1.81/im respectively. Using 0.7" slit-widths (pixel 
scale is 0.18"), the spectral resolution of i? ~ 3600 corre¬ 
sponds to a dispersion of 1.086, 1.303, and 1.627 A/pixel 

IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 


for YJH. The MOSFIRE spectroscopic field of view is 
3' X 6/ and we observed four masks with approximately 
25-35 objects targeted in each to map the galaxy cluster 
(Fig.[T]). All the masks were designed using the MAGMA 
slitmask design software B A slit star was included in 
each mask to monitor the flux throughput and seeing. 
Flat-fields and arcs (Neon, Argon) were taken during af¬ 
ternoon calibrations for each mask. 

For the December 2013 run, we designed three masks 
that were observed with both J and H. The seeing 
FWHM on the first night was ^ 0.7 — 0.9" and improved 
to ^ 0.6" on the second night. Adjusting for the seeing 
conditions to reach approximately the same flux limit in 
all masks, we integrated for a total of 2880 — 3360 sec¬ 
onds and 2880 — 5880 seconds in the J and H bands re¬ 
spectively; individual exposure times for both JH bands 
were 120 sec. All spectroscopy was taken with an ABBA 
dither pattern and an offset of 2.5". Long-slit spec¬ 
troscopy of the standard star Feige 15 was taken on both 
nights to correct for telluric absorption by the atmo¬ 
sphere and to flux calibrate the data. 

For the February 2014 run, we used a single mask to 
obtain Y band spectroscopy of mostly confirmed cluster 
members. Conditions on the first night were poor with 
seeing FWHM of 1.2" (none of these observations are 
usable) but improved to 0.7 — 0.9" for the following 
three nights. The Y mask was observed each night for 
5040 — 6660 sec with individual exposure times of 180 sec. 
The total integration time is 24300 sec, but we note that 
useful time on target is less due to varying conditions on 
the first two nights. Data taken under poorer weather 
conditions were given less weight when combining the 
observations, i.e. the data were weighted by the S/N of 
the continuum. 

To reduce the MOSFIRE spectroscopy, we use the 
publicly available data reduction pipeline (DRP) devel¬ 
oped by the instrument teanB. The DRP provides 
background-subtracted, rectified, and (vacuum) wave¬ 
length calibrated 2D spectra for each slit. Using custom 
IDL routines, we then correct the 2D spectroscopy for 
telluric absorption and flux calibrate using the standard 
star observations. To extract the ID spectra and associ¬ 
ated Icr error spectra, we fit a Gaussian profile along the 
spatial direction to determine the extraction aperture for 
each object; typical extraction apertures are ~ 1". For 
objects that are too faint to fit a spatial profile, we use 
the Gaussian profile ht to the slit star in the same mask. 

The 2D spectra were visually inspected to identify all 
potential sources; for the most part, these detections 
are emission lines, e.g. Ha in the H-band for galax¬ 
ies at z ^ 1.6. Spectroscopic redshifts were determined 
with two-step method: 1) using xcsao in IRAF with an 
emission-line template and, where possible, combining 
JH observations, to measure a redshift and 2) using this 
redshift to ht ID Gaussians to emission lines along the 
spectral direction and propagating errors using the cor¬ 
responding error spectrum. 

For more details on our data reduction, we refer to our 
MOS FIRE analysis of the GOSMOS cluster ([Yuan et al.l 
l2014fl . Especially thorough reviews of MOSFIRE observ¬ 
ing strategies and data reduction are also presented by 

^® www2.keck.hawaii.edu/inst/mosfire/magma.html 

https://github.com/Mosfire-DataReductionPipeline/MosfireDRP 
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iKriek et'all (|2015D and iSteidel et al.l (I2014D ; note that 
the latter uses the public data reduction pipeline (DRP) 
released for MOSFIRE. Both publications provide exten¬ 
sive discussion on the capabilities of the instrument and 
strategies for post-processing. A full description of our 
data reduction for IRC 0218 and the COSMOS galaxy 
cluster at z = 2.1 also will be presented in Nanayakkara 
et al. (in prep). 


2.4. Redshift Catalog 

Each spectroscopic redshift is assigned a quality flag 
Qz where Qz = 3 denotes a robust measurement (multi¬ 
ple spectral lines), Qz = 2 is likely (single spectral line 
with potential secondary line), and Qz = 1 is guess (sin¬ 
gle line and/or no strong spectral features). The com¬ 
bined spectroscopic observations from LRIS and MOS¬ 
FIRE yield 130 unique redshifts. Note that all of the clus¬ 
ter galaxies reported here are confirmed with emission¬ 
lines, e.g. the split [Oii]A 3727A doublet or Ha-|-[Nii] 
pair. Faint continua are detected for only a handful 
of cluster galaxies and are not reliable for measuring 
absorption-line redshifts, i.e. our survey is strongly bi¬ 
ased towards active galaxies with emission lines. In our 
analysis, we use only the 109 unique objects with redshift 
quality flag of Qz = 3 (54 identified by LRIS and 55 by 
MOSFIRE). 

We find the spectroscopic redshifts are in good agree¬ 
ment with the photometric redshifts (Fig. [U top left). 
The (^phof-^spec) distribution for the 99 objects with 
•Zspec> 0 and Qz = 3 has a median difference of [Azjmed = 
—0.006; we use the median to minimize the effect of 
redshift outliers. The corresponding absolute median 
deviation is ^ 3% where the latter is determined us¬ 
ing Az/(1 -I- Zspec) and assuming Zspec=1.62. For the 
49 galaxies at 1.5 < 2 :spec< 1-7, the median difference is 
0.06 and absolute median deviation is also ~ 3%; there is 
only one spectroscopically confirmed member with l^phot- 
^specl ^ 0.3. 

The redshift distributions in Fig. [T] (bottom) demon¬ 
strate how well LRIS and MOSFIRE complement each 
other. LRIS is effective up to z ~ 1.6 because [On] is 
still within the LRIS wavelength range. MOSFIRE is 
extremely effective at z > 1.5 when [Oiil] and Ha move 
into the J and H bands, respectively. 


2.5. NIR Emission-Line Fluxes 

To obtain spectro-photometric data and correct for 
the slit-loss, we flux calibrate the MOSFIRE observa¬ 
tions using a combination of ground-based and Hubble 
Space Telescope imaging. We use the standard star flux 
calibrated data from our custom IDL routines to com¬ 
pare with the da ta from the (publicly ) available UKIDSS 
Data Release 8 (|Ouadri et al.l [201^ . More information 
on the hlter s ystern s used in UKIDSS can be found in 
iHewett et al.l ( 200611 . 

We mask the sky regions in the spectra by assigning 
them weights of zero and use an inverse variance weight¬ 
ing to calculate the total spectroscopic fluxes for each ID 
spectrum. We use the UKIRT WFCAM filter response 
functions to be able to directly compare our results with 
the total broadband photometric fluxes of UKIDSS. 

For each MOSFIRE mask in a given band, we calculate 
the median offset between the spectroscopic magnitudes 


and the broadband photometric magnitudes for all ob¬ 
jects brighter than 23 mag. We use the median offset 
as the correction factor to be applied to all the objects 
in the respective masks to account for slit losses. We 
then calculate the normalized median absolute deviation 
for objects with a broadband magnitude brighter than 
24 mag to determine the uncertainty for the scaling pro¬ 
cess. We find the J-band to have systematically more 
scatter compared to the H-band. We calculate the flux 
uncertainties in the calibrated spectra to be ~20%. A 
more detailed analysis will be presented in Nanayakkara 
et al. (in prep). 

Our analysis uses HaA6563A, H/3A4861A, 

[Nii]A 6583A, and [Oiii]A 5007A line fluxes to mea¬ 
sure redshifts, star formation rates, gas attenuation, and 
gas metallicities of z ^ 1.6 galaxies. Ha is the strongest 
emission line and falls in the H-band; all of the Ha 
detections are measured at > Aa significance for these 
galaxies. We identify 45 galaxies at 1.58 <Zspec< 1-68 
(Fig.|2|) with Ha fluxes of 0.6—16.5 xl0“^^ erg s“^ cm“^. 
The H-band includes the much weaker [Nil] emission 
line; objects with [Nil] measurements that are below 
the detection limit of ~ 0.3 xl0“^^ erg s“^ cm“^ 
and/or detected with statistical significance of < 2a are 
assigned Icr upper limits from the line-fitting. 

For our analysis of star formation and gas metallici¬ 
ties, we use the narrow redshift range of 1.6118 <Zspec< 
1.6348 defined by our kinematic analysis f N3.1l) and ex¬ 
clude the cluster galaxies that are X-ray detected (1) and 
that host broad-line AGN (1; identified by Ha and [Nn]); 
this leaves 26 star-forming cluster galaxies with Ha fluxes 
and a stellar mass range of log(M*/MQ)= 9.2 — 10.9 (Ta¬ 
ble [ij . We consider measurements for both individual 
galaxies as well as for stacked spectra divided into three 
stellar mass bins. 

To stack the spectra, we first divide the 26 star-forming 
members into three stellar mass bins with ranges of 
log(M*/M 0 )= 9.2 — 9.6 (9 members), log(M*/M 0 )= 
9.6 — 10.2 (9 members), and log(M*/M 0 )= 10.2 —10.9 (8 
members). We normalize the individual g alaxies by their 
Ha fl ux to give them equal weight (e.g. ISanders et al.l 
I2015f) . convert the spectra to rest-frame values, combine 
the spectra, and then measure the line fluxes by htting 
Gaussian profiles. Gonfidence intervals for the line fluxes 
are determined by boot-strapping the ID spectra. The 
stacked spectra show even the weaker emission lines H/3 
and [Nil] as well as the [Sii] doublet (Fig. [3]). 


2 . 6 . Stellar Masses & Stellar Attenuation 

To determine stellar masses (M*) and stellar atten¬ 
uation (Ay.sta.r.S B; Table [Bi we use the catalog of the 
UDS field from iQuadri etal.l ([201211 with our spectro¬ 
scopic redshifts. Both M^, and Ay .st.ar.SR are determined 
by using FAST (IKriek et al.ll200g ) to fit Spectral Energy 
Distributi ons (SED s ) to t he multi-wavelength imaging. 
We use a IChabri^ ( 200311 initial stellar mass function, 
constant solar metallicity, and exponentially declining 
star formation rate (r =10 Myr to 10 Gyr). 

We also use FAST to determine the correspond¬ 
ing errors for M* and Ay,star,SB ; the l atter assumes 
Rv=4.05 (starburst attenuation curve; iCalzetti et al.1 
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Fig. 1. — The photometric redshifts from lWilliams et al.l II2009I) compared to spectroscopic redshifts for all objects with measured redshifts 
(top left; the dashed diagonal line denotes parity) and their spatial distribution (top right). Objects with 2:spec quality flag of Qz = 3 are 
shown as open blue triangles (LRIS) and open red squares (MOSFIRE). Note that 30 objects have 2:spec from both spectrographs. The 
bottom panel shows the 2spec histograms for LRIS (blue) and MOSFIRE (red) for objects with Qz = 3. 


l2000fF1 . FAST determines confidence intervals by us¬ 
ing Monte Carlo simulations and modifying the observed 
fluxes by their photometric errors in each simulation. 
Note that there are degeneracies between stellar attenua¬ 
tion, star formation histories, metallicities, and dust laws 
calculated by FAST. These degeneracies can sometimes 
be broken by independently measuring, e.g. nebular at¬ 
tenuation. For a subsample of the galaxies, we measure 
both Ha and H/3 emission to determine the Balmer decre¬ 
ment and thus compare nebular to stellar attenuation 

(B- 

3. ANALYSIS 

3.1. Galaxy Cluster Kinematics 

Our combined LRIS and MOSFIRE spectroscopic sur¬ 
vey confirms 109 unique objects with redshift quality 
flag Qz = 3; of these, nine are stars. We identify 

The starburst (SB) attenuation curve is commonly referred to 
as the Calzetti law and is appropriate for continuum measurements. 
We use “starburst” as requested by D. Calzetti. 


45 galaxies at 1.58 <Zspec< 1-68 and the galaxy clus¬ 
ter is a well-defined peak in the redshift distribution at 
z ~ 1.6 (Fig. [U bottom). To determine the cluster red¬ 
shift and redshift range, we first need to identify the 
members. We remove outliers iteratively using a me¬ 
dian absolute deviation (MAD); the MAD is less sensi¬ 
tive to outliers that can be problematic for dynamically 
young systems such as IRC 0218. For a gaussian distri¬ 
bution, a ss 1.48xMAD and we use 3 ct limits to define 
the redshift range. We employ 5 iterations and find a 
stable solution requires only 2 iterations: IRC 0218 has 
33 members within 1.6118 <Zspec< 1.6348 which corre¬ 
sponds to limits defined by ±3cr (Fig. [5]). 

To determine IRC 0218’s redshift, we bootst rap the 33 
galaxies (15000 realizations: iBeers et al.lHO^ and mea¬ 
sure a median cluster redshift of Zci = 1.6233±0.0003 
and mean cluster redshift of 1.6232 ± 0.0004. The clus¬ 
ter galaxies have a strikingly narrow redshift distribu¬ 
tion and the cluster velocity dispersion (from the boot¬ 
strapped distribution) is correspondingly low: ac\= 254± 
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spec 


X (Mpc) 


Fig. 2.— The galaxy cluster IRC 0218’s redshift range (left panel) is determined by bootstrapping the 45 galaxies at 1.58 <2:spec< 1.68. 
We measure the cluster’s median redshift and velocity dispersion to be z^i = 1.6233di0.0003 and cr^i= 254 ± 50 km s“^; the latter is 
shown as a Gaussian curve. The cluster redshift range is defined by di3c7ci (left panel, dotted vertical lines) which isolates 33 emission-line 
members (filled histograms). The spatial distribution of the confirmed cluster members is shown in the right panel (filled circles); galaxies 
outside the cluster redshift range (open histogram) are shown as open circles. At z = 1.62, the projected angular scale is 8.471 kpc per 
arcsecond; the large dotted circles correspond to projected (proper) radii at 2 = 1.62 of 0.5 and 1 Mpc. 


TABLE 1 

Properties of Hq-detected IRC 0218 Members^ 


Keck ID 

a(2000) 

5(2000) 

I mag 

^spec 

log(M*/MQ) 

-^V,star,SB 


fH/3 

'b 

f[Nii]b 

SFR(H«3tar)‘= 

36395 

34.61690 

-5.20299 

24.70 

1.613 

9-4itg:J| 

Q 2+0.3 

U.Z_o 2 

1.8 

± 

0.1 

0.0 

± 

0.6 

0.4 

± 

o 

w 

2.2 

36849 

34.56660 

-5.19900 

25.00 

1.624 

9.60tni 

0 1 

3.4 

± 

0.2 

0.9 

± 

0.4 

0.2 

± 

0.2 

3.6 

37269 

34.56870 

-5.19590 

24.72 

1.624 

10.39tg-5f 

n 1+^-^ 

5.1 

± 

0.5 

0.4 

± 

0.5 

0.8 

± 

0.3 

5.3 

37522 

34.60300 

-5.19369 

24.95 

1.622 

W.20t°om 

-^•-*^-0.5 

2.1 

± 

0.2 

0.4 

± 

0.3 

0.7 

± 

0.1 

8.1 

38080 

34.61000 

-5.18868 

24.15 

1.625 

10.34t°°“ 

O.+a? 

10.3 

± 

0.2 

2.6 

± 

0.5 

1.6 

± 

0.2 

11.0 

38455 

34.60940 

-5.18625 

24.02 

1.623 

10.87t°;°| 

0.8lO;i 

4.6 

± 

0.2 

1.2 

± 

0.3 

0.8 

± 

0.2 

12.2 

39150 

34.59990 

-5.17829 

25.50 

1.624 

y.zz_o 06 

U.^_0.2 

2.3 

± 

0.1 

0.3 

± 

0.5 

0.4 

± 

0.1 

2.8 

39463 

34.59290 

-5.17625 

24.23 

1.622 

Q c-t+0.05 

y-o<_o.io 


00 

± 

0.2 

1.7 

± 

0.6 

0.7 

± 

0.1 

9.2 

39771 

34.59280 

-5.17402 

24.36 

1.623 

9.811^6 

U.O_0.2 

4.9 

± 

0.3 

1.3 

± 

0.4 

1.6 

± 

0.2 

10.1 

39989 

34.58750 

-5.17219 

25.28 

1.623 

10.29t°;J“ 

^•^-0.2 

3.6 

± 

0.3 

0.0 

± 

0.3 

2.5 

± 

0.3 

4.3 

40243 

34.61680 

-5.16958 

24.78 

1.622 

9-6ilS:SI 

U.^_0.2 

2.4 

± 

0.2 

0.3 

± 

0.5 

0.4 

± 

0.1 

2.8 

40382 

34.61750 

-5.16884 

25.18 

1.623 

10.05t°;?| 

-] 4+0-2 
-*-•^-0.2 

5.3 

± 

0.2 

1.1 

± 

0.3 

1.1 

± 

0.2 

30.8 

40568 

34.58070 

-5.16690 

24.79 

1.626 

y-^y_o.05 

n f)+0-3 

6.7 

± 

0.3 

1.4 

± 

0.2 

0.8 

± 

0.1 

6.2 

40729 

34.58480 

-5.16558 

24.76 

1.627 


fl 1 

2.0 

± 

0.4 

0.2 

± 

0.4 

0.0 

± 

0.1 

2.1 

40731 

34.59280 

-5.16529 

24.84 

1.625 

y.zo_o.i3 

n f)+0-3 

2.7 

± 

0.1 

2.2 

± 

0.7 

0.5 

± 

0.2 

2.5 

41189 

34.62220 

-5.16206 

23.77 

1.623 

y- <o_0.05 

0-2la2 

9.0 

± 

0.5 

2.8 

± 

0.3 

1.4 

± 

0.2 

10.9 

41297 

34.60190 

-5.16097 

24.44 

1.622 

9.93t“;“ 

O.slO;^ 

5.3 

± 

0.3 

0.5 

± 

0.4 

0.8 

± 

0.1 

14.0 

41548 

34.59790 

-5.15948 

24.48 

1.624 

I0.50t°;p 

1 q+0-3 

7.3 

± 

0.3 

0.1 

± 

0.5 

2.0 

± 

O 

CO 

25.1 

41956 

34.59960 

-5.15566 

23.75 

1.621 

10.25t°;“® 

1 fl'^^’^ 
-*^•'^-0.2 

14.6 

± 

0.4 

2.6 

± 

00 

o 

4.3 

± 

0.2 

50.2 

43178 

34.61020 

-5.14441 

24.66 

1.621 

I0.47t°;« 

T Q + 0.6 

-*-•^-0.3 

5.4 

± 

0.3 

1.1 

± 

0.5 

2.0 

± 

0.2 

24.3 

44102 

34.56320 

-5.13659 

24.62 

1.621 

10.69t°?o 

fl 

13.8 

± 

o 

CO 

1.0 

± 

0.6 

4.1 

± 

0.2 

21.5 

44403 

34.57880 

-5.13336 

24.88 

1.622 

9.36l°;JJ 

fl 2^^'^ 
^•^-0.2 

4.9 

± 

0.3 

2.0 

± 

0.6 

0.3 

± 

0.1 

5.9 

44459 

34.62260 

-5.13383 

23.93 

1.624 

10.121012 

n 4+0-3 

7.0 

± 

0.3 

0.2 

± 

0.5 

0.5 

± 

0.3 

11.0 

44587 

34.56660 

-5.13347 

24.12 

1.621 

y-y^_o.i4 

n c + 0.2 
^•^-0.1 

4.9 

± 

0.3 

2.2 

± 

0.5 

0.5 

± 

0.1 

8.8 

45959 

34.56070 

-5.12105 

25.05 

1.620 

Q OO + 0-22 

y.4»_Q Q6 

l-elas 

2.0 

± 

0.2 

0.3 

± 

0.4 

0.0 

± 

0.1 

15.2 

46922 

34.61150 

-5.11372 

23.21 

1.630 

iU.iO_Q Q7 

0.4l0;0 

14.9 

± 

0.4 

4.2 

± 

0.4 

5.7 

± 

0.3 

23.6 


^ IRC 0218 members (1.6118 <2:spec< 1.6348; see ^3.111 that have Ha emission as measured with MOSFIRE. 

^ Observed fluxes in units of 10“^^ erg s“^ cm“^. Note th at in several cases, fH/d and f[Nii] are upper limits considering the S/N. 
^ Ha fluxes are corrected for dust using Ay,star (Eq-[Zl s®® ^3.211 and then converted to star formation rates using the relation from 
IHao et ^ Il2011lj . 




















7 



03 



6400 6500 6600 6700 6800 

a[A] 

Fig. 3.— We divide the 26 HQ;A6563A-detected cluster galaxies into three separate mass bins and stack their spectra; for clarity, the 
stacked spectra are smoothed and offset in flux. The stellar mass range is log(M*/MQ)= 9.2 — 10.9 (TableO and the median stellar mass 
for each bin is included in the figure; from low to high mass bin, there are 9, 9, and 8 members. The individual spectra are first normalized 
by the Ha flux and converted to rest-frame values, and then the median is taken. At Zd = 1.6233, the spectral range corresponding 
to MOSFIRE includes H^A4861A and the [Oiii]AA4959, 5007A pair (top; J-band), and Ha, [Nii]A 6583A, and the [Sii]AA6717, 673lApair 
(bottom; H-band). 
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E(B-V) 


star,SB 


Fig. 4.— Left: Distribution of Ha vs. H/3 line fluxes ( xlO erg s ^ cm as measured with MOSFIRE where the diagonal line 
corresponds to Ha/H/3= 2.86 (Balmer decrement for Case B recombination). Our requirement of H/3> 0.5 xl0~^^ erg s“^ cm“^ measured 
at > 3(7 confidence identifies 12 cluster galaxies (filled circles); included are the Icr errors determined from propagating the error spectra. 
The attenuation for the the 3 cluster galaxies with Ha/H/3< 2.86 is likely negligible and so we set their Balmer decrement to 0. Right: 
The color excess measured using the Balmer decrement compared to the color excess from the stellar continuum measured with FAST; we 
include only the 9 cluster galaxies with Ha/H^> 2.86 and use the CCM89 reddening law to calculate the nebular E(B-V). For reference, 
we include the ratio of stellar t o nebular E(B-V) measured for local starburst galaxies where the color excess for stars is less (/ = 0.44, 
dashed line dCalzetti et al.ll200C)h . The cluster galaxi es at ^ ^ 1.6 follow this general trend, but there is signicant scatter that mirrors recent 
results based on larger samples of field galaxies fsee lReddv et a 012013). 


50 km s“^, i.e. more like a group than a fully devel¬ 
oped galaxy cluster (see 1 14.1I1 . Figure [2] (right) shows 
that most (20/33) of the cluster members are within 
Rproj< 1 Mpc of the BCG. If we consider only the mem¬ 
bers at Rproj< 0.5 Mpc, the cluster redshift and veloc¬ 
ity dispersion are virtually identical. The combination 
of high spatial density and low velocity dispersion are 
very c onducive to galax y-galaxy merging (|Rudnick et al.l 
120121 : iLotz et iini2013D . i.e. frequency and duration of 
galaxy-galaxy interations is higher than in a field enyi- 
ronment. 

IRC 0218’s redshift and yelocity dispersion do not 
change significantly if we limit our analysis to either 
the MOSFIRE or LRIS obseryations only. If we repeat 
our analysis using only the MOSFIRE redshifts, there 
are 30 members: the median cluster redshift remain the 
same and the velocity dispersion decreases slightly to 
CctMOS = 220 ± 43 km s“^. Using only the LRIS red- 
shifts, there are 17 members: the median cluster redshift 
remains the same and the velocity dispersion increases 
slightly to fTci,LRis = 305 ±118 km s“^. 

3.2. Balmer Decrement and Attenuation 

For a subset of MOSFIRE-detected galaxies, we are 
able to measure both H/3 and Ha and thus directly deter¬ 
mine their internal extinction via the Balmer decrement. 
Ha is by far the stronger line (Fig. [3]) where Ha/H/3=2.86 
is the intrinsic line flux ratio for Case B recombination 
with electron temperatu re Tp, = 10^ K and electron den¬ 
sity of Ue = 100 cm“^ (lOsterbrockI 119890 ■ In addition 
to being weaker, at 0 ~ 1.6 the H/3 line falls in a region 
of telluric absorption (Aobs ~ 1.264//m) which makes it 
even more challenging to measure. 

Only 12 cluster galaxies meet our selection criteria of 
H/3> 0.5 xlO”^^ erg s“^ cm“^ measured at > 3cr sig¬ 


nificance (Fig. m left). There are 9 cluster galaxies with 
Ha/H/3> 2.86 and three members with lesser values that 
are within ~ 1 — 2(T of this ratio (Fig. HI left). The latter 
are likely to have negligible attenuation and so we set 
their attenuation to zero. 

The total attenuation curve is parameterized by 


p - 

^^-E(B-V) 

(1) 

k{X) = 

^ ’ E(B-V) 

(2) 


where A(A) is the total magnitude of the extinction at 
wavelength A and E(B-V) is the color excess as measured 
by nebular lines or the stellar continuum. For galaxies 
with a measured Balmer decrement, the color excess is 
determined directly from the observed line fluxes with 

E(B - V)hii = [ k{m - MHa)]-' 2.51ogio 

,_. (3) 

Using the ICardelli et al.l (CCM; 119891) attenuation 
curve for the diffuse interstellar medium where Ry = 3.1, 
/c(Ha) = 2.53, and fc(H/3) = 3.61, the color excess is then 

E(B - V)hii = 2.33 log^o • (4) 

Figure 0] (right) compares the color excess of the ion¬ 
ized gas E(B-V)hii to that of the stellar population E(B- 
V)star; the latter is determined by fitting to the (UV) 
stellar continuum wi th FAST (solar meta llicity) and as¬ 
suming Rv,sb=4.05 (jCalzetti et al.ll2000l) . For the error 
in E(B-V)hii, we add the Ha and H/3 line flux errors in 











































9 


quadrature. For the cluster galaxies, the nebular E(B- 
V)hii is larger than that of the stellar component E(B- 
V)star,SB- This is consistent with results at z ~ 0 and 
z ~ 1.4 that support a ratio of E(B-V)=t.a,r to E(B-V)hii 
of 0.44 (ICalzetti et al.l[2000l : IPrice et al.fl2014D . However, 
we note the large scatter in this ratio that mirrors results 
from recent studies of field galaxies at z ~ 2 (iReddv et al.l 
1201511 . Although it is tempting to measure this ratio us¬ 
ing the cluster galaxies, our limited sample does not allow 
us to confirm nor exclude the ratio of 0.44. 

Figure [S] (left) shows the cluster galaxies that satisfy 
our rigorous H/3 criteria and compares their nebular at¬ 
tenuation to stellar mass. The members do not define 
a correlation between A(Ha;)HTT and in contrast to 
results at z ^ 0 .1 (iGarn fc Best! 1201011 and z ~ 1.4 
(jPrice et al.ll2014ll . nor is any trend visible when we stack 
the spectra into three mass bins. 

The cluster galaxies define a clearer trend between stel¬ 
lar attenuation Ay, star,SB as measured by FAST (Fig. [5l 
right) and M*. The least-squares fit to the 20 member^ 
is 

Av,star.SB = 0.48 log[MfyMQ] - 4.36 (5) 

where ctrms = 0.34. A fit using only the 9 members 
with a measured Balmer decrement is identical. The re¬ 
lation between Ay^star.SB and M* see ms to be shallowe r 
than measured for the field at z ~ 1.4 (jPrice et al.ll2014l l. 
but there is large scatter. 

How then can we correct our Ha fluxes for attenua¬ 
tion given that for most of the Ha-detected galaxies, we 
only have the stellar attenuation Ay^star.SB measured by 
FAST? For simplicity, we first correct the Ha line fluxes 
using the nebular attenuation curve from ICardelli et al.l 
(fT^ with Rv=3.1 

A(Ha)Hii = 2.53 x E(B - V)h„. (6) 

We then use the observed stellar to nebul ar attenuation 
ratio of E(B-V)star= 0.44xE(B-V)hii (jCalzetti et al.l 
l2000fl and the stellar attenuation Ay^star.SB as mea¬ 
sured by FAST for all of th e galaxies where E(B- 
V)star=Av,star,SB/4.05 (see also lSteidel et al.| |2Ml) 

A(Ha)Hii = 5.75xE(B-V),,,, (7) 

We use Eq. [3 to correct all of the Ha fluxes for attenu¬ 
ation an d determine corresp onding star formation rates. 
Recentlv lReddv et al.l (1201^ measure Rv=2.505 for field 
galaxies at z ~ 2, a value that is lower than the Ry for 
the starburst and diffuse ISM attenuation curves. How¬ 
ever, we cannot discriminate between these models given 
our limited cluster sample. 

3.3. Ha Star Formation Rates 

To determine Ha star formation rates, we correct for 
attenuation using E(B-V)star (Eq.[7]) and convert the Ha 
fluxes to s tar fo rmation rates using the relation from 
iHao et a].l (l^fTTill : 

log[SFR(Hastar)] = log[L(Hastar)] - 41.27 (8) 

We exclude X-ray and broad-line AGN from the fit and also 
member 45959 whose SED fit indicates a very dusty low-mass 
galaxy (see Table [T}. 


This relatio n assumes a Kroupa IMF (0.1 — 100 Mq; 
iKroupal 120011) . but the relation for a Chabrier IMF is 
virtually identical (difference of 0.05). The uncorrected 
Ha fluxes of the 26 star-forming members range from 

I. 8 — 14.9 xlO”^"^ erg s“^ cm“^. Using E(B-V)star to 
correct for attenuation, the corresponding Ha star for¬ 
mation rates are 2 — 50 Mq yr“^ (Tabled]). We have 
excluded X-ray detected members and broad-line AGN. 

The cluster field was imaged with Spitzer MIPS as part 
of the legacy UKIDSS Ultra Deep Survey (SpUDS PI: 

J. Dunlop) and the catalog includes all sources detected 
with S/N > 5 which corresponds to a flux of ^ 40/iJy. 
The 24/im fluxes are c onverted to tota l infrar ed luminosi¬ 
ties (L/r) using the iGharv fc Elbad (120011) templates, 
and S FRs are ca l culate d from Lin using the prescrip¬ 
tion of iKennicuttI (|1998[ ). adjusted to the Ghabrier IMF. 
We are able to measure 24/rm fluxes for the six cluster 
members with the highest Ha (> 20 Mq yr“^). 

Even with correcting for dust attenuation, the Hastar 
SFRs are lower compared to the 24^m SFRs: on av¬ 
erage, the 24^m SFRs for the six detected members 
(42 — 155 Mq yr“^) are a factor of ~ 3 larger. The 
advantage of using Ha is that we are significantly more 
sensitive, i.e. we can measure SFRs to 2 Mq yr“^ com¬ 
pared to ^ 40 Mq yr“^ with the 24/im observations, and 
thus we are not limited to only the mos t strongly star¬ 
formin g dust-obscured cluster galaxies. iPapovich et al.l 
(120071) also showed that f or galaxies at z ^ 2, SFRs based 
on iGharv fc ElbiQ ()2001[ ) templates fit only to 24/im are 
larger by ~ 2 — 5 compared to SFRs based on a com¬ 
bination 24/im, 70/im, and 160/im. This is intriguingly 
consistent with the offset we measure when comparing 
our Hastar SFRs to the 24/im SFRs. 

We confirm the original measurement by iTran et al.l 
(|20ir)t) of an elevated star formation rate within the clus¬ 
ter core: the integrated Hastar SFR per square Mpc 
at Rproj< 0.25 Mpc is about six times higher than at 
Rproj< 2 Mpc (Fig. ini Tabled]). The five star-forming 
members at Rproj< 0.25 Mpc have stellar masses of 
log(M*/MQ)= 9.3 — 10.3. The cluster galaxies also show 
increasing SFR with stellar mass (Fig. [ 7 ]). We stress that 
the integrated Hastar SFRs are lower limits because 1) 
for the six members with both Hastar and 24/im SFRs, 
the 24/im values are on average ~ 3 times higher; and 2) 
the limited spatial coverage of the NIR spectroscopy com¬ 
bined with the increasing contamination by field galaxies 
means we identify fewer Ha-emitting members at larger 
Rproj- In M.41 we compare in detail the cluster Hastar 
SFRs to that measured for the field at z 2. 

3.4. Gas-Phase Metallicities 

Following seve ral recent analyse s on gas metallici¬ 
ties in the field (iKulas et al.l 120131 iSteidel et al.l 12014 
I Sanders et ahlPM^ . we measure the oxygen abundance 
for the cluster galaxies u sing the N2 method calibrated 
bv IPettini fc Pag^ (120041) : 

12 -I- log (0 /H)j^ 2 = 8-90 + 0-57 x log ([Nii]/Ha) (9) 

We use only N2 to measure gas metallicities due to 
our limited ability to measure H/3 and [Oiil] for indi¬ 
vidual galaxies. We exclude X-r ay detected members 
and broad-line AGN; we also follow fSanders et all (|2015D 
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Fig. 5. — The nebular (left) and stellar (right) attenuation versus stellar mass for individual cluster galaxies with measured Balmer 
decrement (filled circles) and Ay, star,SB only (open circles). Shown for comparison are the stacked values for the nebular (left) and stellar 
(right) attenuation of the field at 2 1.4 (open stars: [Price et aT1l2014l) . The IRC 0218 galaxies do not show any trend between nebular 

attenuation A(Hq;)hii and stellar mass; this is also true when stacking the cluster spectra. In contrast, the cluster galaxies do increase in 
stellar attenuation Ay, star,SB with M* as confirmed with a least-squares fit (diagonal line), but the cluster has a shallower slope than in 
the field (see also lReddv et ani2015ll . 




Rp.„i (Mpc) (Mpc) 


Fig. 6. — The integrated Hostar star formation rate per unit area (the total Mf:; yr~^ per square Mpc) is hig hest in the cluster core 
(Rproj< 0.25 Mpc; left) and decreases with increasing cluster radius fe.g. I'ltan et al.ll2010l : IRrodwin et al.l[2013ll . This is driven by the 
existence of star-forming members in the core, but the average Hostar SFR per galaxy is actually lowest in the cluster core and plateaus 
at Rproj^ 1 Mpc (right; Table [2j. The propagated errors due to the uncertainty in the Ha line fluxes are smaller than the points 


(±0.6 - 1.4 Mq yr-l). 


and require that log([Nii]/Ha)< —0.3 to remove poten¬ 
tial AGN which removes one member. From our sample 
of Ha-detected cluster galaxies, we measure N2 for the 
20 members where the [Nil] emission line is detected at 
(S/N)> 2; note that [Nil] is much weaker than Ha. For 
five additional members, we place Icr upper limits on N2. 
We also measure gas metallicities for the binned spectra 
(Fig.©. 

The cluster galaxies show increasing metallicity with 
mass (Figs. |8] & ©, a trend that is consistent with 


the mass-metallicity relati on (MZR) measured fo r star¬ 
forming galaxies at z 0 (jMoustakas et al.ll2011h . The 
cluster galaxies lie below the local MZR, i.e. the clus¬ 
ter galaxies at z ~ 1.6 have lower gas metallicities at 
a given stellar mass. The offset in the cluster MZR is 
similar to t hat observed for field star-for ming galaxies 
at z 1.6 (IZahid et al.||2 014[) at z 2.3 (ISteidel et al.l 
12014 lenders et al.ll2015lb 


4. RESULTS 
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TABLE 2 

Hcustar Star Formation Rates vs. Cluster Radius 


Rproj 

Ngal 

S(Ha SFR)^ 

S(Ho SFR) per galaxy 

<SSFR> <= 

S(H 

a SFR) per Area 

250 kpc 

5 

28 Mq yr“^ 

5.6 Mq yr“^ 

0.8 Gyr-i 

144 

Mq yr“l Mpc“^ 

500 kpc 

8 

62 Mq yr“^ 

7.8 Mq yr-l 

0.9 Gyr-i 

79 

Mq yr“l Mpc“2 

1 Mpc 

16 

197 Mq yr-1 

12.3 Mq yr ^ 

0.8 Gyr-i 

63 

Mq yr“i Mpc“^ 

< 2 Mpc 

26 

324 Mq yr“^ 

12.5 Mq yr“^ 

0.9 Gyr-i 

26 

Mq yr“^ Mpc“^ 


^ The total enclosed Ha star formation rate corrected for attenuation using E(B-V)star (see ^3.211 . 
X-ray detected members and broad-line AGN are excluded. Given spectroscopic incompleteness at 
larger radii, we consider these to be lower limits on IRC 0218’s total SFR. 

^ The integrated Hagtar SFR divided by the number of Ha-detected galaxies within Rproj ; quiescent 
members are excluded. 

^ The specific star formation rate determined by summing the Hagtar SFR and M* for the Ha- 
detected members within Rproj; quiescent members are excluded. 




Fig. 7.— The 26 cluster galaxies tend to have lower Hagtar star formation rates (left) and specific SFRs (right) at a given stellar mass 
compared to their field counterparts. Left: The Hagtar SFRs for the cluster galax ies (filled circles) are corrected for extinction using 
E(B-V)gtar and typical errors in Hagtar SFRs and M* are ~ 0.1 — 0.2 dex (see ^ 13.21 Table^. Included are lines of constant specific star 
formation rates (SSFR; Mq yr“^ per M*; dotted lines) and the horizontal dashed line is our low er SFR limit. For co mparison, we include 
the binned points measured by MOSDEF for Ha-detected field galaxies at 2 : 2.3 (open triangles: ISanders et al.ll2015ll and the star-forming 

field galaxies at 1.5 < 2 : < 2.0 from ZFOURGE (long-dashed curve; Tomczak et al, submitted; UVH-MIR SFRs). The ZFOURGE errors 
represent uncertainty in the mean and the stellar mass bins are 0.25 dex. Right: The higher mass cluster galaxies tend to have lower SSFRs 
compared to their field counterparts in both ZFOURGE (long-dashed curve) and MOSDEF; a least-squares fit to the IRC 0218 members 
with SFRs> 2 Mq yr“^ (solid line; la range shown as thin lines) measures a slope that is steeper than in ZFOURGE. The dashed diagonal 
line shows the SSFR limit corresponding to our SFR limit of 2 Mq yr“^. 


4.1. Is IRC 0218 A Galaxy Cluster? 

IRC 0218 is one of only a handful of galaxy clusters at 
z > 1.5 with a measured cluster velocity dispersion ad. 
The challenge lies in obtaining medium-resolution near- 
infrared spectroscopy for a sufficient number of members 
to measure accurate redshifts and map the cluster’s kine¬ 
matic structure (see 13.ip . Although NIR spectroscopy 
is strongly biased towards emission-line galaxies due to 
the faintness of the stellar continuum, this may be ad¬ 
vantageous giyen__theinnreased_star_Jormntion at this 
epoch (e.g. iTran et alJ[20ir)HBrodwin et alJl201,‘^ . Note 
that while cluster studies at ^ ^ 2 using HST/WFC3 
grism observations can confirm both active and passive 
members via, e.g. th e D4000 break (|Gobat et alJ 120131 : 
iNewman et alJ [20141 . the grisms’ redshift accuracy of 
(1 + x) X 1000 km s“^ leads to rather uncertain measure¬ 


ments of the cluster velocity dispersion, e.g. iMei et ahl 

dMl). 

Using Ha to id entify 15 emission-line members, 
iBavliss et al.1 (I2014D measure a very large velocity dis¬ 
persion of 1500 ± 520 km s“^ in a galaxy cluster at 
z = 1.48. In comparison, the velocity dispersions of 
~ 750 — 850 km s ~^ measured for two p roto-clusters 
at z = 2.2 and 2.5 (|Shimakawa et all 1201411 are more in 
line with those of lower red shift clusters, e.g. EDisCS 
(jMilvang- Jensen et al.ll200^ . The smallest velocity dis¬ 
persion observed thus far at z ^ 2 is for the COSMO S 
cluster with ad = 552 ± 52 km s“^ (jYuan et al.ll20l4l : 
this cluster also has the most confirmed members (> 50) 
as identified by our ZFIRE survey. 

Relative to these studies, IRC 02I8’s velocity disper¬ 
sion of ad= 254± 50 km s“^ (33 emission-line members; 
Fig. [2]) is by far the smallest. This low velocity dispersion 
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Fig. 8.— Stellar mass versus gas metallicity for IRC 0218 members as measured by [Nii]/Ha (filled red circles); we fit a least-squares to 
these galaxies (red line). These cluster galaxies are also divided by SFR[Hostar]=10 Mq yr“^ into low (small filled red circles; 9) and high 
(large filled red circles; 11) star formation rates. Red open circles (5) with arrows are memb ers with Icr upper limits on [Nil]. Included 
for comparison are: the MZR of field galaxies at ^ ^ 1.6 in COSMOS (blue long-dashed curve: IZahid et al.ll2014|): field galaxies at 2 : ^ 2.3 
from MOSDEF (open green tria ngles:TSanders et al.ll2015lh the MZR relation from SDSS (dashed black curve: iMoustaka s et al.||2011|): the 
linear fit to 2 2.3 galaxies from ISteidel et al.l (dot-dashed cyan line: 1201^ : and solar abundance (dotted line at 8.69: F^phm^^i^ni2009l) . 

The MZR relation for the cluster IRC 0218 (red line) is clearly offset from local galaxies (dashed curve) and agrees with the Steidel fit to 
galaxies at 2 ~ 2.3 (dot-dashed cyan line). 


is more in the rang;e of KO-lo-xy groups than a g alaxy clus- 
ter (e.g. lGerke et al.l[2005l : iTanaka et alJl2013i) . However, 
this velocity dispersion is consistent with X-ray results: 
the cluster velocity dispersion corresponding to the total 
cluster mass of M 200 = 7-7 ± 3.8 x 10^^ Mq estimated 
from the weak X-ray si gnal detected with XMM-Newton 
(2.3 ct level detection; iPierre et al.l I2012H is CTci,xray 
360 ± 90 km s“^, i.e. the dynamical velocity dispersion 
is consistent with the X-ray derived estimate at ~ l.Str. 
These values also are con sistent with pr e dictio ns from nu¬ 
merical simulations, e.g. lEvrard et al.l (j2008[ l show that 
a halo with X-ray mass of M 200 7.7 x 10^^ Mq has 

an estimated velocity dispersion of 400 km s“^. Note 
that the cluster members are spatially concentrated with 
20 of the 33 members within Rproj= 1 Mpc of the BCG 
(Fig. [2l right). 

Using the Millenium simulation with a semi-analytic 


galaxy formation model, iShattow et al.l (I2013D select 50 
central galaxies in halos with masses greater than 7 x 
10^^ Mq and track them from z = 6 to z = 0; note 
that this halo mass is essentially IRC 0218’s dynamical 
mass as measured by both kinematics and X-ray. The 
study finds that by z = 0, most of the halos have masses 
> 2 X 10^"^ Mq, i.e. they are lik ely to be clusters if 
not the most massive ones (see also lChiang et al.ll2013H . 
Thus IRC 0218 will evolve into a g alaxy cluster that can 
be as massive as the Virgo cluster (|Schindler et al.l[T999H 
and possibly as massive as 10^® Mq. 

4.2. Nebular vs. Stellar Attenuation 

E(B-V)hii and E(B-V)star can be used to determine 
dust properties, i.e. E(B-V)hii probes sightlines to re¬ 
gions of ionized gas that are heated by short-lived 
OB stars while E(B-V)star is a measure of the (UV) 
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Fig. 9. — The same as in Fig. |8] but for the stacked spectra where th e cluster galaxies a r e divided into three m ass b ins (large r ed fille d 
circles; see Fig. [S)!. We show the field data binned in stellar mass from IZahid et aP II2014I 1. [Sanders et al.l II2015I V and lErb et al.l Il2006al l. 
The IRC 0218 members are remarkably consisten t with the field MZR relation at 2 ~ 2 as well as the MZR measured for the COSMOS 
cluster (filled orange circles: [Kacorzak et al.ll201Rll . 


dust attenuation of the stars and integrated over longer 
timescales (> 100 Myr). For local starbursts, E(B-V)star 
and E(B-V)hii are correlated such that the color excess 
of the stellar co ntinuum is about hal f that of the ionized 
gas (/ = 0.44: iCalzetti et aI1l2000H . i.e. the attenua¬ 
tion as measured by the nebular lines is larger. It is 
unclear if this relation holds at z > 1 with studies mea- 


suring ratios varving from ^ 0.26 — 1 (lErb et al.ll2006bl: 

Yoshi 

tawa et al.l 

20101: IKashino et al.l 120131: IZahid et al.1 

ml' 

Price et al. 

Iml IReddv et a.lJl201,5n. 

We 

are able to directly measure the color excess E(B- 


V)hii for 9 cluster galaxies where Hq;/H/ 3> 2.86 and 
H^> 0.5 xlO“^^ erg s“^ cm“^ is detected at > Scr 
(Fig. m right). Both E(B-V)hii and E(B-V)star span 
ranges (0 — 0.8) that are comparable to values mea- 
sure d at z ^ 1 — 2 in both field and c luster stud¬ 
ies (jZeimann et al.l 120131 [Price et al.l l2014[l . The ra¬ 
tio of nebular to stellar attenuatio n is also consistent 
with the measured val ues of ^ 0.44 (ICalzetti et al.ll2000l : 
lYoshikawa et al.ll2010f) . albeit with considerable scatter. 


There is no clear trend between the nebular attenua¬ 
tion A(Ha)Hii and stellar mass (Fig. [5]), nor is there any 
indication that nebular attenuatio n depends on environ¬ 
ment as suggested by past studies (iKovama et al.l[^13bl : 
IZeimann et al .1120131) . In contrast, the stellar attenuation 
Ay, star,SB (starburst, Rv=4.05) does increase with stel¬ 
lar mass (Fig. [5]) but with a slop e that is shallowe r than 
observed in the field at z 1.5 (jPrice et al.ll2014ll . Our 
results suggest that the attenuation measured by the in¬ 
tegrated stellar light may be less affected by, e.g. patchy 
dust distribution; howe ver, our sample siz e is limited. 
We refer the reader to iReddv et al.l l|2015ll for a more 
extensive study of dust based on a significantly larger 
sample of (field) galaxies at z ~ 2. 

4.3. Star Formation in the Cluster Core 

With our sensitive spectroscopic survey of IRC 0218, 
we measure extinction-corrected Hogtar SFRs for 26 in¬ 
dividual members with ^2 — 50 Mq yr“^; AGN have 
been removed. Note that an advantage of spectroscopy 
over using narrow-band imaging to derive Hogtar SFRs 
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is that we can remove AGN using [Nil]/Ha as well as 
those detected using X-ray. With the 26 cluster galax¬ 
ies, the lower limit on IRC 0218’s integrated Hogtar star 
formation rate is ~ 325 Mq yr“^ (Table |2|). The higher 
integrated Hogtar SFR in the cluster core (Fig. jS]) sup¬ 
ports our earlier result bas ed on 24/rm observations of 
IRC 0218 (iTran et al.ll2'M^ . 

A handful of studies have now also obtained Ha spec¬ 
troscopy of cluster galaxies at z ~ 1.5 using ground- 
based telescopes and the WFC3 grism in Hubble Space 
Telescope. Line sensitivities vary, but the consensus is 
that most galaxy clusters in this redshift range have 
integrated SFRs within Rproj< 0.5 Mpc of up to ^ 
200 Mq yr“^ and SF Rs of individual members can be as 
high as 100 Mq yr“^ (iHavashi et al.l[ 20 TillZeimann et al.l 
1201311 . With an Hagtar SFR of ~ 60 Mq yr“^ integrated 
within Rproj< 0.5 Mpc (Table [2]), IRC 0218 is consistent 
with these measurements. 

Our analysis uses the integrated Hagtar SFRs within 
different Rproj (Fig. [U Table |2|) to enable direct compar¬ 
ison to the literature. Note that the integrated Hagtar 
SFR within Rproj is different from the fraction of star¬ 
forming members. Because the overall galaxy density is 
higher in clusters, it is possible to have both a large inte¬ 
grated SFR and a n increase in quiesce nt gal axies relative 
to the field, e.g. iQuadri et al.l (|2012l l and iSantos et al.l 
([201411 . For example, the Hagtar-SFR per square Mpc is 
highest in the core (Rproj< 0.25 Mpc), but the average 
SFR per galaxy is only about half that of galaxies at 
Rproj < 2 Mpc (Fig. mi . 

Our results confirm that the integrated Hagtar SFRs 
in the cores of galaxy clusters at z > 1.5 are signif- 
ican tly higher than in massive galaxy clu s ters at z < 
1.4 (Kodama et ^12004 iFinn et ^ 120051 : I Bauer et al.l 
1201 lib The astute reader will notice that our Hagtar 
SFRs are l ower than the 24/rm values published in 
ITran et ^ ([201011 . For the six IRC 0218 members with 
both Hagtar and 24/rm SFRs, the average ratio is ^ 
2.9± 1.1, i.e. the 24/rm SFRs are ^ 3 times higher. How- 
ever, it is likely that t he 24/xm SFRs are over-estimated: 
iPapovich et al.l ( 200711 showed that for galaxies at z ~ 2, 
SFRs based on Chary &: Elbad ([200111 templates fitted 
only to 24/im are ^2 — 5 higher compared to SFRs de¬ 
termined with 24/rm, 70/im, and 160/im. 

We note that if we use the 24/im measurements for 
the six Ha-detected members, the the number of cluster 
galaxies with Ha SFRs> 30 Mq yr“^ increases from two 
to six. Because we are also constrained by the spatial 
sampling of the NIR spectroscopy, we consider the Hagtar 
SFRs in Table|5]to be lower limits. However, this result is 
tentative and a larger cluster sample at z 2 is needed 
to better quantify the relation between Ha and 24/rm, 
especially for SFR> 20 Mq yr“^. 

4.4. Star Formation Rate versus Stellar Mass 

The 26 Ha-detected cluster galaxies have stellar masses 
of ~ 10® — 10^^ Mq and show a trend of increasing SFR 
with stellar mass that is similar to the field (Fig. [71 
left), i.e. a “main sequence” of star-for ming members 
([Daddi et al.ll2007t iKovama et al.l 1^13bll . Most of the 
lower-mass cluster galaxies (M*< 10^® M©) '^iU more 
than double their stellar masses given the elapsed cos¬ 
mic time since z = 1.62 if they maintain their current 


Specific SFRs (SS FR; see also Table ip . Our results 
are consistent with iZeimann et ‘all ([201311 who use WFC3 
grism measurements of Ha to show that cluster galaxies 
at 1 < z < 1.5 follow a SFR-M* relation. 

To determine if the SFR-M* relation depends on 
environment, we us e the Four-Star Gala xy Evolution 
(ZFOURGE) survey ([Tomczak et al.ll^0l4 Straatman et 
ah, submitted) to select star-forming field galaxies at 
1.5 < z < 2.0 with SFRs measured by combining UV 
and mid-IR fluxes Tomczak et ah, submitted). Given our 
Hagtar detection limit, we apply a SFR detection thresh¬ 
old of 2 M 0 yr"^ to both the IRC 0218 and ZFOURGE 
samples. We compare Hagtar SERs and the specific star 
formation rate (SSER defined to be SER/M*) to stellar 
mass in Eig. jT] We also compare to recent results from 
MOSDE E based on Ha meas urements of field galaxies at 
z rv. 2.3 ([Sanders et al.l 1201^ . 

The IRC 0218 members have systematically lower 
SFRs compared to the field. However, the difference 
can be due to using two diffe rent methods of measurin g 
SFRs (Hagtar vs. UV-I-MIR: iKennicutt fc Evansll20l3l . 
Instead, we consider the slope to test for any environ¬ 
mental difference. By fitting a least-squares to IRC 0218 
and bootstrapping the sample ( 10 ^ realizations), we mea¬ 
sure the SSFR-M* slope to be —0.51 ± 0.30. For the 
ZFOURGE field galaxies, we measure a shallower slope 
of —0.23 ± 0.07. The difference in the SSFR-M* slope 
between IRC 0218 and ZFOURGE suggests that envi¬ 
ronment plays a role, but our results are not statistically 
significant. 

An alternate comparison is to iSanders et al.l (120151 1 
who measure Ha SERs for the field at z 2.3 (Fig. [71 
The IRC 0218 members are clearly offset to lower SFRs 
compared to the Ha-detected field galaxies. One possible 
explanation for the lower SFRs in IRC 0218 is the elapsed 
cosmic time between z ~ 2.3 and z ^ 1.6 {At 1.1 Gyr). 
However, the offset in the SFR-M* relation for MOSDEF 
field galaxies at z ~ 1.6 at z ~ 2.3 (iKriek et al.l 1201 511 is 
not as significant as observed in IRC 0218. We conclude 
that our observations indicate a dependence on environ¬ 
ment, i.e. lower SFRs in the cluster especially for the 
higher mass members. 

4.5. Gas Phase Metallicity-Mass Relation (MZR) 

There is considerable and confusing disagreement over 
the observed relation between gas phase metallicity and 
stellar mass (MZR) at z 2. While several studies of 
field and protocluster ga laxies find that metallicities in¬ 


crease with stellar mass (iZahid et al.ll2014l: ISteidel et al. 

2014 

IShimakawa et al.ll2015ll. others do not (Stott et al. 

2013 

iKulas et al.l2013l). Using the N2 method calibrated 


bv lPettini fc Pagell ([20041 1. we measure the oxygen abun¬ 


dance in cluster galaxies and compare to their stellar 
masses (Figs.| 8 ]&|^. We find that the IRC 0218 members 
are offset to lower metalliciti es relative to star-for ming 
galaxies in the local universe ([Tremonti et al.ir2004f) . 

The cluster MZR is re markably similar to that of sev¬ 
eral f i eld surveys at z ^ 2 ([Erb_et_ajj2^r)06a ySteidel et al.l 
12014 lenders et al.ll2015HZahid et al.ll2014l l: this is par¬ 
ticularly eyident when comparing the stacked spectra 
(Fig. |S]). To quantify the trend of increasing metallic¬ 
ity with stellar mass, we fit a linear least-squares to the 
20 members with measured [Nil]/Ha (we exclude the fiye 
members with only upper limits on [Nil]): 
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12 + log(0/H)N2 = 8.49 + 0.13[log(M*/Me)-10.0] (10) 

where ctrms = 0.10. The MZR in IRC 0218 is com¬ 
parable to the relation measured bv iSteidel et al.l (l2014f l 
for field galaxies at z ~ 2.3 where the slope and normal¬ 
ization are (0.2,8.41). The consistency holds across the 
full range in mass [9 <log(M*/M0)< 11]. For these 20 
members, the median difference bet ween measured [Nii ] 
and that predicted using the fit from iZahid et al.l (|2014f ) 
for field galaxies at z ^ 1.6 is only 0.04. 

Whether there is an environmental imprint on the 

MZR at z ~ _ 2 is yet to be clearly established. 

lEllison et al.l (|2009H detect a small offset of ~ 0.04 dex to 
higher metallicities in overdense environments at z ~ 0 
(see also lCooner et al.ll2008l : [Peng fc Maiolinoll2014[ l. and 
recent studies at z 2 claim to find a measurable offset 
betw e en their cluster and fie ld samples (iValentino et al.l 
12014 iShimakawa et all l2015f) . However, we find no evi¬ 
dence of an environmental dependence for the MZR when 
comparing IRC 0218 to the field at z ~ 1.6. 

What is striking is how the cluster MZR does not 
evolve at 1.6 < z < 2.1 (over ~ 1 Gyr): IRC 0218’s MZR 
is virtu ally identical to the C OSMOS cluster at z = 2.1 
fFig.l^ lKacprzak et al.ll2015li . The similarity of the clus¬ 
ter MZRs matches the essentially constant g as metallic¬ 
ities of leased field galaxies at 1.5 < z < 2.5 (lYuan et al.l 
l2013fl . Our results are also in line with th e slow evo¬ 
lution of the MZR predicted by simulations (iDave et al.l 
l2011al) , but we note that these same simulations also pre¬ 
dict an environmental imprint, e.g. higher metallicities 
at a given stellar mass, that we do not observe in either 
galaxy cluster. 

Our results show that the cluster MZR and its evolu¬ 
tion is indistinguishable from that of the field at 1.6 < 
z < 2.3 which is consistent with predictio ns from recent 
simulations bv lTavlor fc Kobavashil (I2015D . We conclude 
that the somewhat contradictory results observed thus 
far for the MZR at z ^ 2 are due to a combination of 
different line flux sensitivities and sample sizes of the 
various studies. By expanding our sensitive survey of 
cluster gal axies at z ^ 2 and comparing to results from 
MOSDEF (|Kriek et al.ll2015h . we hope to resolve these 
disagreements. 

4.6. Star Formation Rates and the MZR 

Studies at low redshift find that at a given stel¬ 
lar mass, the gas metallicity is lower for galaxies 
with higher star formation rates (iMannucci et al.l 120101 : 
[Andrews &: Martinil I2013D . The Fundamental Mass- 
Metallicity relation (FMR) quantifies this dependence of 
the MZR on the SFR and can be explained by a SFR de¬ 
pendence on the ratio of self-enrichment versus inflow of 
pristine gas. Whether the FMR holds at z ~ 2 is unclear 
with recent stu dies finding no evidence of an MZR depen - 
dence on SFR (iWuyts et al.ll20l4 iSanders et al.l[20T^ . 
Evidence for an environmental depend ence for the EMR 
at z 2.3 is also tentative at best (|Shimakawa et al.l 

nnn. 

To test for a dependence on SFR, we divide the 
IRC 0218 members into approximately equal sized low 
and high SFR bins (SFR[Hastar] = 10 Mq yr“^; Fig.|S]); 
we exclude the five members that have only upper lim¬ 
its on their gas metallicities because their true values 


could be considerably lower. At high stellar masses 
[log(M*/M0)> 10], the high SF members dominate with 
a metallicity range of ~ 0.2 dex. At low stellar masses, 
the low SF members dominate with a larger metallicity 
range of ~ 0.4 dex. 

We fit a least-squares to the 11 high SFR members 
and the 9 low SFR members with measured [Nil] and 
measure slopes with errors of 0.05 ± 0.08 and 0.10 ± 0.11 
respectively. The MZR for both low and high SFR bins 
are consistent (within Icr) with that of the c luster sample 
as a w hole as well as the MZR measured bv iSteidel et al.l 
(|2014[) . We conclude that the IRC 0218 members show 
no evidence of an MZR dependence on the SFR, but we 
acknowledge that our cluster sample is limited. 

5. EVIDENCE OF ENVIRONMENTAL DEPENDENCE? 

With our NIR spectroscopy, we measure Ha and [N ii] 
fluxes that, combined with our SED fits from FAST, pro¬ 
vides extinction-corrected Ha star formation rates, gas 
phase metallicities from [Nil]/Ha, and stellar masses. 
The primary evidence we find for any environmen¬ 
tal dependence is in the cluster star formation rates. 
IRC 0218’s integrated Hastar SFR per unit area is about 
a factor of three higher at Rproj< 0.25 Mpc compared to 
Rpr,^ 1 Mpc (from 144 to 26 M0 yr“^ Mpc“^; Ta¬ 
ble [2]). This is driven by the higher number density 
of star-forming galaxies in the cluster core relative to 
studies of cluster co res at lower redshifts (|Tran et al.l 
120071 : iMei et al.l[20n9[ i. However, the average Hastar SFR 
per galaxy is actually lower in the core compared to 
the values at Rnrni= 0.5 and 1 Mpc (Table [U see also 
iBrodwin et al.l[201^ . 

At a given stellar mass, the IRC 0218 members have 
Hastar SFRs that are lower compared to the field at 
z 2.3 (Fig. [7]). The slope of the SFR-M* relation 
also seems to depend on environment such that the mas¬ 
sive cluster galaxies have lower SFRs compared to their 
field counterparts. The quenching of SF in the mas¬ 
sive cluster galaxies at z ^ 1.6 means these members 
will pas sively age to ma tch their cluster counterparts at 
z < 1 (IMei et al.l 1200911 . The continued stellar growth 
for lower mass cluster galaxies is also consistent with the 
more extended star fo rmation histories measu red in clus¬ 
ter galaxies at z < 1 (iTran et al.ll2005al 1200711 . 

Like field galaxies at z ~ 1.5, the cluster galaxies show 
an increasing amount of attenuation with stellar mass, 
i.e. more massive galaxies are dustier. This trends hold 
primarily for the stellar attenuation (Ay,starve) as mea¬ 
sured by the SED fitting to the multi-wavelength pho¬ 
tometry (Eig. O. However, the large scatter in both 
in IRC 0218 and field studies (lYoshikawa et al.l 120101 : 
iPrice et'mi2014ji is sobering, particular when compar¬ 
ing the nebular attenuation A(Ha)Hii to M*. It is very 
likely that the empirical relations used to determine at¬ 
tenuation and thus obtain extinction-corrected SERs at 
z ^ 2 have un comfortably large u ncertainties and depend 
on SERs ('e.g. lReddv et al.ll2015ll . A better characteriza¬ 
tion of dust attenuation laws for difference spectral types 
at z > 1 using essent ially low resolution spectroscopy 
(jKriek fc Conrovl[2013l l holds promise in addressing this 
issue. 

The gas phase metallicity versus M* (MZR) for 
IRC 0218 members is nearly identical to the field at 
z ~ 1.6 (Eigs. [8] & [9l). The cluster MZR does not de- 
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pend on SFR which is also consis t ent w i th field results 
(IWuyts et alJ 12014 iSteidel et alJ 12014 I Sanders et alJ 
l2015li . Compa ring to our results for the COSMOS clus¬ 
ter at z = 2.1 (iKacorzak et al.l[201!Tf) shows no evolution 
in the MZR in the ^ 1 Gyr. Our ZFIRE cluster results 
contradict recent claims at z ^ 2 of a measurable differ¬ 
ence in the MZR with environment (|Kulas et al.l 1201^ 
iShimakawa et al] 1201511 . We note that the environmen- 
tal imprint on the MZR is small even using SPSS at z ^ 0 

(iCooper et al.ll200MlEllison et al.l[2009l:lPeng fc Maiolinol 

I20l4) . thus if there is an environmental dependence for 
the MZR at z ^ 2, a considerably larger sample of cluster 
galaxies is needed to detect it. 

Alternatively, calibration issues with using [Nil]/Ha 
to measure gas phase metallicities at z ~ 2 may mask 
any environmenta l impr int on the MZR. For example, 
iKewlev fc EllisonI (j2008[ l showed that the choice of line 
ratios used to measure the gas-phase metallicity has a 
significant effect on the shape of the MZR. By obtain- 
ing the full suite of e mission lines for the BPT diagram 
([Baldwin et al.l HoMI) to separate star-forming galaxies 
from those hosting ACN, we can compare to metallici- 
ties using alternativ e line ratios such as H/3/[Oiii] (e.g. 
iShaolev et al.l[2C)l^ . We can also characterize the ion¬ 
ization conditions. If there is a relationship between 
SSFR and ionization parameter, this would have inter¬ 
esting implications for the properties of star clusters at 
high redshift, e.g. the spatial distribution and m asses 
of stars embedded within them (|Kewlev et al.ll201^ . By 
expanding our survey of cluster galaxies at z ~ 2, we will 
explore how gas metallicities and ionization parameters 
are linked to star formation, galaxy (stellar) mass, and 
environment. 


6. CONCLUSIONS 

We survey the galaxy cluster IRC 0218 using the 
multi-object spectrographs LRIS (optical) and MOS- 
FIRE (near-IR) on Keck I. We obtain 130 unique red- 
shifts based mostly on emission lines, e.g. [On] and Ha, 
and identify 45 galaxies with 1.58 <Zspec< 1.8 (Figs.[T]& 
HI). The IRC 0218 members define a narrow redshift dis¬ 
tribution with 33 galaxies at 1.6118 <Zspec< 1.6348. The 
corresponding cluster redshift and velocity dispersion are 
Zci = 1.6233 ± 0.0003 and CTci = 254 ± 50 km s“^. 

While IRC 0218’s velocity dispersion is the lowest mea¬ 
sured thus far for a cluster at z ~ 2, it is consistent 
with the a inferred from the cluster mass estimate of 
M200 = 7.7 ± 3.8 X 10^^ M (7i from the weak X-ray de¬ 
tection (jPierre et al.ll201^ . Simulations also show that 
most systems with this velocity dispersion and redshift 
will evolve into a cluster by z ~ 0. IRC 0218 is spatially 
concentrated with 20 spectroscopically confirmed mem¬ 
bers within Rproj< 1 Mpc of the BCG. The cluster’s 
high spatial concentration and low velo city dispersion 
are conducive to galax y-galaxy merging ([Rudnick et al.l 

IWaiFntz emil20Tli . 

With MOSFIRE, we reach near-IR line flux sensitiv¬ 
ities of 0.3 xl0“^^ erg s“^ cm“^ that, combined with 
our SED fits to multi-wavelength photometry, provide 
extinction-corrected Ha star formation rates, gas phase 
metallicities using [Nil]/Ha, and stellar masses. Adopt¬ 
ing a Hastar SFR limit of 2 Mq yr“^ at Zd = 1.6233 
and removing ACN, we confirm our earlier result based 
on 24/xm observations of elevated star formation in 


IRC 0218’s core (jTran et al.ll2010D : the integrated SFR 
per unit area at Rproj< 0.25 Mpc is about a factor of 
three larger than at Rproj< 1 Mpc (Fig|6l Table [2]). We 
consider the Hastar SFRs to be lower limits due to sparse 
sampling and because for the handful of objects where 
we also can measure the Balmer decrement and/or 24^m 
flux, the individual SFRs increase by ~ 2 — 3. 

In our analysis, the strongest evidence for any en¬ 
vironmental dependence at z ^ 1.6 is in the cluster 
star formation rates. The high Hastar SFR in the core 
(R proj< 0.25 Mpc; Fig. H]) is driven by the concentra¬ 
tion of star-forming members while the average SFR per 
galaxy in the core is actually half that of galaxies at 
Rproj"^ 1 Mpc (Table [2]). The massive cluster galaxies 
also tend to have lower SFRs compared to their field 
counterparts (Fig. [7|). This mass-dependent quenching 
of SF is needed for these cluster galaxies to match the 
stellar ages of their cluster coun terparts at z < 1 (e.g. 
iTran et al.l[200^ IMei et al.l[200ll . 

Like the field, the IRC 0218 members show a trend 
of increasing attenuation with stellar mass as measured 
by the SED fits (Ay star, sb), be. more massive galax¬ 
ies are dustier (Fig. [^. However, this is not the case 
when using the nebular attenuation as measured by the 
Balmer decrement (A(Ha)Hii)- The significant scatter 
in the nebular attenuation suggests larger problems with 
correctly measuring dust attenuat ion at z 2 for b oth 
cluster and field galaxies (see also lReddv et al.l[2015fl . 

The gas phase metallicity versus stellar mass (MZR) 
for the cluster galaxies is virtually identical to the field 
at z ^ 1.6, and both are offset to lower metallicities com¬ 
pared to galaxies at z ^ 0 (Figs. pfcip ). Comparing to 
the fit determined bv iZahid et al.l ( 20141 for field galaxies 
at z 1.6, the median difference between measured and 
predicted [Nil] is only 0.04. Fitting a least squares to the 
IRC 0218 members confirms that th eir MZR is the same 
as for the field as measured by, e.g. ISteidel et al.l (|2014D 
at z ~ 2.3. The cluster MZR does not depend on star 
formation rate which is also consistent with recent re¬ 
sults on field galaxies (jWuvts et al.l 1201 J : iSanders et al.l 

[Ml). 


Particularly striking i s that IRC 0218 and the COS¬ 
MOS cluster at z = 2.1 ([Kacorzak et al.ll2015f) have the 
same MZR even though the elapsed time is ~ 1 Gyr. 
Both clusters’ MZRs are also well-matched to the field 
as measured by several surveys at z ~ 2. We do not see 
any evidence of an environmental impr int on the MZR at 
z ^ 2 as indicated by p revious studies (|Kulas et al.ll2013l : 
IShimakawa et al]l2015f) . We attribute this difference to 
our combination of a larger cluster sample with a more 
sensitive line flux limit. 

In summary, environmental effects are likely to be sub¬ 
tle at z ~ 2 and require a larger survey of cluster galax¬ 
ies to fully characterize. With ZFIRE, we are building 
such a cluster sample and also expanding to include addi¬ 
tional line diagnostics needed to test, e.g. the reliability 
of [Nil]/Ha for measuring gas phase metallicity. ZFIRE 
complements ongoing NIR sp ectroscopic surveys of field 
galaxies at z ^ 2, e.g. KBSS ([Steidel et al.ll2014) . MOS- 
DEF (iKriek et al.ll201^ . and KMOS^° ([Wisnioski et al.l 
Ml, and will enable us to test for evolution in spatially 
resolved quantities as a function of environment. 
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